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Abstract 

The Casimir energy, free energy and Casimir force are evaluated, 
at arbitrary finite temperature, for a dilute dielectric ball with uniform 
velocity of light inside the ball and in the surrounding medium. In 
particular, we investigate the classical limit at high temperature. The 
Casimir force found is repulsive, as in previous calculations. 



1 Introduction 

Calculations of the Casimir energy for spherically symmetric boundary con- 
ditions have been the issue of numerous papers, since the first calculation 
by Boyer [|l| (see also [0, H, H, ||, ||). Although a general theory describing 
the fluctuations of the electromagnetic field at finite temperature has been 
present for a long time []^, |, ||, |1T|, the extension of the results for the 



dielectric ball from zero temperature to finite temperatures received less at- 
tention. The case of a conducting spherical boundary (at zero and finite 
temperature) was studied by Balian and Duplantier p[, and will serve here 



to verify the validity of the results we obtain. Brevik and Clausen [0 , used 
the Debye asymptotic expansion to calculate the Casimir force on the dilute 
dielectric ball, assuming uniform velocity of light inside and outside the ball 
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(including, however, frequency dependence of the permeabihty and permit- 
tivity), to leading order in the so-called dilute approximation (the nature 
of this approximation is explained following Eq.(|TB|)). It turns out that the 
usual Debye approximation is most useful for zero temperature calculations, 
but involves complications at finite temperatures. For example, it was shown 



T3| that use of these asymptotics may lead to problematic results, such as 
that the Casimir energy of a dielectric ball is independent of temperature. 

In this Letter we calculate the thermodynamic features of the Casimir ef- 
fect for the dilute dielectric ball at finite temperature (to leading order in the 
dilute approximation). We obtain explicit simple expressions for the Casimir 
energy, free energy and for the Casimir force at arbitrary temperature. In 
particular, we find that (within the approximation made) the Casimir force 
is repulsive for any radius at any temperature (see Eq.(^)). At high tem- 
perature its leading behavior is proportional to jg^, where a is the radius 
of the ball. Higher order corrections to the dilute approximation decay like 
powers of ^, and consequently the leading behavior gives the Casimir 
force quite accurately at high temperature (see section 5). 

The starting point of the present calculation is an explicit expression 
for the Casimir energy density (to leading order in the dilute approxima- 



tion). This expression was derived in using simple properties of the 
Green's function of the Helmholtz equation. The Casimir energy at finite 
temperature, Ec{T), is then obtained by summing that explicit expression 
for the energy over the Matsubara frequencies. Then, using the usual ther- 
modynamic relation between the Casimir free energy Fc{T) and the Casimir 
energy Ec{T), we obtain an explicit expression for Fc{T) valid for all tem- 
peratures. We verify that EciT) and Fc{T) coincide at T = 0, as they 
should (and the Casimir entropy vanishes at T = 0). Furthermore, at high 
temperatures, Fc{T) tends to a term which is proportional to — Tlog(aT) 
(where a is the radius of the ball). This behavior is consistent with the re- 
sults of and is also consistent with the recent general discussion of the 
high-temperature classical limit in |TH]. 



2 Mode summation at finite temperatures - 
general considerations 

The eigenfrequencies of the electromagnetic field, subjected to some bound- 
ary conditions, are solutions of a set of characteristic equations 

A,^(c.)=0 A; = 1,2,3... (1) 
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labeled by an index k. (The eigenfrequency spectrum is the union of solutions 
to all the equations in (0).) The set of parameters c describes the constraints, 
such as the radius of the sphere or the distance between the conducting plates. 
The functions A^(ci;) are usually expressions involving solutions of the scalar 
Helmholtz equation subjected to appropriate boundary conditions. In the 
following we consider the case of spherical symmetry, where the functions 
A^(co') involve combinations of Bessel functions (see Eq. ( |l^ - |T6|) in sec. 3). 

Let us first consider the system in its ground state at zero temperature. 
The zero temperature Casimir energy is of course the difference between the 
formally divergent sum over all allowed modes of the constrained system and 
the divergent sum over the modes of the free electromagnetic field, 

n n 

This expression is still divergent and needs regularization. Thus, we introduce 
a UV-frequency cut-off and define 

When the boundaries are pushed to infinity, the system ceases to be con- 
strained and becomes free. This is described by letting the parameters c in 
d^) tend to an appropriate value, say, infinity. It is convenient to represent 
the sum in the last expression as a contour integral 

Ec{u:N) = Y.i ^Fk{z;c)dz (4) 

k "^^iV 

where r^r is shown on Fig.^, and 

F,{z;c)= hm -L ± (log A'^{z) - log A^z)) , (5) 

c'->/ree ZlTl dZ 

is the (regulated) resolvent corresponding to index /c in (|lD, and "/ree" de- 
notes the set of parameters describing the limiting case where the boundaries 
are taken to infinity. In the case of spherical geometry this is achieved by 
taking the radius to infinity. The Casimir energy is then defined as 

Ec = lim Ec{ujn)- (6) 

To carry over these considerations to finite temperature, we recall that the 
average energy per mode of the electromagnetic field is ^ coth(/3a;/2) and 
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Figure 1: A rectangular contour in the complex plane, 
that the free energy per mode is log(2 sinh(/3co'/2)), where P = -j^. Thus, 
Ec{T-un) = Y.I ^-coi\i{^)F,{z-c)dz, (7) 

k ^'^N 

is the (cut off) Casimir energy at temperature T, and 

Fc{T; cuN) = J2(f ^ log(2 smh{Pz/2))Fkiz; c)dz, 

k -^^N 



(8) 



is the (cut off) Casimir free energy at temperature T. Upon rotating the 
integration contour in (|^) to wrap around the imaginary axis, the Casimir 
energy can be expressed as a sum over Matsubara modes, namely 



c 



EE(I^' 



n=—oo k 

This concludes our general considerations. 
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(9) 
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3 Casimir energy 



In this section we study the temperature dependence of the Casimir energy of 
a dielectric ball, having permittivity and permeability e, fi which is embedded 
in a medium with permittivity and permeability e', fi'. A further simplifying 
assumption which we make throughout this Letter, is that the velocity of 
light be continuous across the boundary, namely 

e/i = e'/i'. (10) 

This condition is often desirable as it causes the frequency equations to sim- 
plify and some divergences to cancel out (see |jl6[ and references therein). It 



turns out that in this case the energy depends on the dielectric data only 
through the parameter defined by 



which is symmetric under ei ^82- Thus, the Casimir energy is invariant 
under interchanging the inner and outer media. 

It has been shown [|], |l6l that under the condition ([To|), and after ro- 
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fating the integration contour in (^ to the imaginary axis, that the difference 
in zero point energy of two balls of radii a and h may be written as 



Ec{a) - 

1=1 



1 

Ec{h) = Y,{1 + -) / duuFiiu- a, 6), (12) 
where the resolvent (^) is 

&) = -^^ ( log [1 - e>^M^\)\ - log [1 - em^\)\ ) . (is) 

Here 

\i{a\uj\) = {siei{a\u\))\ (14) 
where s;, e/ are the modified Bessel functions (i/ = / + |) 

si{x) = ixji{ix) = y-^e 2 J^[ix) (15) 



ei{x) = ixh] [ix) = '^y ^ H^[ix) (16) 
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The case of ^ finite but small corresponds to the inner ball being nearly 
identical in its properties to the surrounding medium. This is referred to in 
the literature as the limit of dilute media. In this limit we can expand (|TB|) 
in powers of This expansion is further justified by the fact that on the 
imaginary uj axis is small, thus the expansion can be pushed to ~ 1 
(for = 1, of course, one of the media has infinite conductivity and is thus a 
perfect conductor). In this Letter we will content ourselves with the leading 
term in this expansion. Thus, to order we have from (|13]) 



FM, ^) = £^ UMuj\) - \]{h\uj\) ] , (17) 



In one of us was able to carry the sum over I in (|T2D to first order in 
in closed form, without using Debye approximations, as is usually done in 
the literature. The result was 

1=1 

Thus, the relative resolvent, to order is 

F^^\z-a,h) = [e~^''\'\{l + Aa\z\+Aa^\z\^)-e-^''\'\{l + Ah\z\+Ah''\z\^)\ 

(19) 

(here z is a pure imaginary frequency.) 

Using this simple explicit expression, we can calculate the energy in the 
dilute limit at an arbitrary temperature. In order to calculate the zero point 
energy with respect to the vacuum we take the limit 6 — oo at the outset. 
Note that the x independent term in (|18|) cancels between the two terms 
in the difference (plOf), resulting in an expression which decays exponentially 
when 1^1 — »• oo. Thus, no further regularization will be needed when we 
integrate over u). Using (^ and (p^Of ), the Casimir energy at temperature T 
is 

Ec{a, T) = ^ En=-oo e-S"^^"!"! (1 + 8aT7r|n| + IQa^T^n^n^) = (20) 



yielding 



1 — -U -' " O". 1 2 _|_ 



Ec{a, T) = ^ + 2(e8JT„i) + (gBa.T_i)^ + (21) 



(e8a,rT_l)3 
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This result is exact (to order for all temperatures. It is immediate to 
check that as T goes to zero we retain the correct zero temperature result 



BcKO) = |^. (22) 



Note that at high temperatures, T » we have 
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Ec{a,T) = ^ + 0{e-'^-'), (23) 

which is consistent with Eq.(8.39) in for the free energy, to leading or- 
der. Namely, the Casimir energy at high temperatures is independent of 
the radius of the ball. That Ec(T,a) becomes independent of the radius 



at high temperature is expected on general grounds (see [jT5[ for a detailed 
discussion): At high enough temperatures we expect classical equipartition 
to hold. Since there is obviously a one to one correspondence between the 
states of the system with radius a and the states of any other system with a 
different radius b, the difference between the zero point energies of these two 
systems must vanish at high temperatures. Indeed, for the case studied here, 
the difference (to order C,^) vanishes exponentially fast with temperature. 



4 Free energy and force 

Having calculated the energy we may now use the expression (|2T|) to infer 
the Casimir free energy of the system. From the general relation 

we have F = ^ J E{(3)d(3 + ^ (where C is an integration constant to be 
determined), or, equivalent ly, 

F = -T y ^^dT + CT . (25) 

To evaluate the Casimir free enrgy at low temperatures we first expand Ec (T) 
around T = 

Ec(a, T) = — ^ + ^ + ^ + 0(T^) 26 

' ' 32a7r ^ 45 945 ^ v ; v ; 
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(Note that the first temperature dependent contribution to the energy is of 
order T^.) Thus, from (p5|), the Casimir free energy is 

^) - + 135^ 472^ + ) • 

To fix the constant C, we use the third law of thermodynamics, namely, that 
the Casimir entropy 

S,(a,T) = ^2(-2ll^Ik(<R ,27) 

vanishes at T = 0. It is straightforward to check that this implies C = 0. 
Therefore, at low temperatures, the Casimir free energy is: 

^) - 32^ 135^ 472^ + ^ ' ^''^ 

This result coincides with the two scattering approximation for the con- 
ducting sphere (up to multiplication by ^^). 
To proceed beyond low temperatures we rewrite ( pS] ) as 

F = -T r + E(0) (29) 

./n 



One can easily check, using the power series (26), that the integral on the 
right side converges, and that indeed F satisfies the condition (pll). To show 
that there is no additional term of the form CT in (^9]), we note that this ex- 
pression can also be verified for any particular mode uj of the electromagnetic 
field. Writing, 

OJ . r % coth ^ - ^ oj 
Tlog(2sinh— ) = -T / ^ (30) 



and then summing over the modes gives precisely (^91). 

In order to evaluate (12^) it is convenient to split Ec into 



Ec{a,T) = EciiT) + Ec2(T) (31) 

where 

EciiT) = ^ T- H 32) 
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and 



^-(T) = 4 + ^-5^ (33) 



and then split the integral in (P9D accordingly. The contribution associated 
with Eci can be integrated in a straightforward manner, whereas the contri- 
bution from Ec2 requires more careful treatment. Thus, 

^ dt ^ f^.^ , , ^ , ..dt 



Fc = -T (Edit) - EciiO))-, - T / {Ec2{t) - ^C2(0))-^ + 



(34) 

The term associated with Eci is now straightforwardly integrated to yield 



3 5T avrT^ (ST+SaTrT^ 



+ ^ + 



32a7r 16 (gSa^T _ 1)2 g (gSa^T _ 

We are now left with the second term, which can be written as 



(35) 



F2{T) = -T {Ec2{t)--^)-- = -^ 0-7 + 



32a7r't2 2 \2 t - 1 J t 

(36) 

The last integral is, of course, well defined at T = 0, but diverges logarith- 
mically as T — s> 00. To simplify this expression we observe that 

°° ^1 1 1 \ e"^* 1 1 

1 — : dt = log T(z) — z log z ^ z — log z log 27r 

2te*-iyt ^ 2^ 2^ 

(37) 

(eq. 8.341.1 in |^), where 2; > can be thought of as a UV-regulator, which 
is set to zero in the end (it is understood that SairTz « 1 throughout the 
calculations.). Thus, writing 



l-SanT ^ 1 \ 1 



Urn + dt 

z^oj^ \2 t e'-lj t 

and using the expansions 

\ogr{z) = -\ogz--^z + 0{z'^) (39) 
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(6.1.33 in [|2T|], 7 = 0.577... is Euler's constant) and 

-dt = -7 - log(8a7rT2) + OiSairTz), (40) 



SawT t 

(eq. 8.214.1 in [^), we find after taking the hmit z — > 0^, that 

{\-\ + ^) \dt = (41) 

2 log(4aT) + I + — JganT (e'-l)i^^ ' 

in which the last integral is well defined for all T. Combining this result with 
(|35|), we obtain from ( p^ the following expression for the Casimir free energy 
at arbitrary temperature: 

Fc{a, T) = -^ [log(aT) + log4 + 7 - |] + (42) 

2 JsoTrT (e'-l)i '[^sTPrZrf 8{e8"'rl'-l) 

Note that log 4 + 7 — | = 0.714, and thus, the linear and Tlog(aT) terms 
in ( ^2] ) agree with the two-scattering approximation (when ^ = 1). All 
other terms in (^) decay exponentially as T ^ 00 and not as powers as 
is seen in eq. 8.39 0]. This is due to the fact that the density of states in 
the approximation decays exponentially with frequency and not as 1/a;^, 
which was the basis of the high temperature approximation used in [0] (see 
eq. (6.12) therein) . This suggests that at sufficiently high temperatures, 
we will have to include higher order corrections in to see the power law 
decay. However, as we shall see in the next section, there are no further 
corrections to the T and Tlog(aT) terms which depend on the radius a. 
Thus, the approximation gives a reliable estimate of the Casimir force at 
high temperatures. Another interesting feature of the last expression is that 
Fc{T) < for a large, while Fc{T) > for small values of a. Thus, for any 
temperature T there is an intermediate radius b{T) such that Fc = for this 
radius (it follows that the free energy may be also viewed as relative to the 
energy of the ball with this zero-energy radius.) 

From (H2]) we calculate the Casimir force to order 



^cia,T) = -^^lFcia,T) = (43) 
re ^ re ^ r2^2g8a.T ^ 27rT3e2(e8'^-^ + l)e8«-^ 



16a% 8a%(e8'^-^ - 1) a?{e^^^T_iY a(e8"-^ - 1)3 

We see that the force (^) is positive for all values of a and T. Moreover it 
behaves as -,^^^3^ at high temperatures. Another quantity of interest, which 
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is the Casimir entropy (pT]), behaves as ^(1.714 + log(aT)) at high temper- 
atures. 

5 Relation to direct mode summation 

In this section we use a different approach based on evaluating (H) directly. 
This analysis was carried out for high and low temperatures in the case of 
the conducting sphere by Balian and Duplantier . To make our discussion 
well defined we consider the difference in free energy between two balls of 
radii a and h. From we deduce that the relative Casimir free energy is 

"1 * /*00 /O 

Fc(a)-Fc(6) = 5^(/ + -)^ / FK^;a,6)log(2zsin(^))rf.;. (44) 
1=1 ^ •^-'^ 

It is convenient to decompose the logarithm as 

log(22sin(|^)) = (45) 

log (2 |sin(f ) I) + .TT E:=o ^ - ^) - El=o ^ (-^ - 

(where, as usual, ^^^q ^i^ans that the n = term is counted with half- 
weight). Note that Fi{u;a,b) is an odd function of frequencies. Thus, dis- 
carding the even part of (^), we are left with 

Fc{a) - Fc{b) = f Er=iil + I) EZi Fii^; b)du + (46) 

The first term on the right side is a decaying function of temperature (it 
decays as |;). The linear and logarithmic contributions in T come only from 
the second term. 

We now expand (PB]) in powers of and work to leading order in 
Substituting F^'^\uj;a,b) from ([T9|) into (|^), we can show that the leading 
term in ( ^Bj ) is consistent with the expression (^) we obtained in the previous 
section. (In the following we integrate without using the b dependent terms, 
wherever they are not necessary for convergence.) We start with the easily 
integrable expressions: 

I EZi (4ae-^- + Aa^u;e-'ndu: = | EZi ^^^'^^'^^ (47) 

Saiv SaTT 

/32(-l+eT)2 
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Note that we summed starting with n = 1, since the n = term is indepen- 
dent of the radius and thus cancels out (when taking the difference). The 
less trivial terms can be calculated as follows 

j2 



^ n=0 — 

pi J>° poo pb p2 rb °° i 

2/3 ^ J2n. / 2/5 A ^ u 



(48) 



"2^ i^-2^°^y+y^^^ 



oo 



-du - / — — rdu 

1 + e") 78^M(-l + e") 

Thus the difference between the free energy of two balls is: 



Fc{a,T) - Fc{h,T) = 



- log(^) + Mi^-naT) - MiMhT) 
2 b 



(49) 



M{x) = \ \ ^ + I /, ^,,, du (50) 



where 

xe"^ + 5(— 1 + e^'' 

4(-l + e-)' ' J. u(-l + e")' 

Clearly (^9]) and (|^) are consistent. 

Finally, we make the important observation that terms of order and 
higher in the expansion of ( ^6] ) do not contribute to Fc{a,T) terms that are 
proportional to T or Tlog(aT) with coefficients that depend on a. To see 
this, note from (|T3|) and from ( ^Bf ) that contribution to terms proportional 
to T or T log(aT) will come from expression such as 

^,ET=iil + -2)rF!'-\co;a,b)duj= (51) 



(The other terms, namely, terms coming from the first sum in (^), obviously 
decay as powers in ^.) For any m > 1 each term in the last sum can be 
integrated to yield: 



(52) 



since for all m > 1 this sum is convergent at = (note that \i{0) = 2rn)- 
Thus, the T and Tlog(aT) terms which depend on a in the exact Fc{a,T) 



are accurately accounted for by the order expression (^). Consequently, 
I) gives the Casimir force quite accurately at high T, up to corrections in 



1 

T- 
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6 Conclusions 



The temperature dependence of the Casimir energy, free energy and force 
were studied in the dilute dielectric ball approximation. We showed that the 
free energy, and thus the force, can be obtained from the Casimir energy, 
which is an easier quantity to calculate. This was done for an arbitrary 
temperature, and agrees with the two-scattering results for a conducting shell 
obtained in for the high and low temperature limits. We confirmed that 
in the high temperature limit the free energy behaves as — Tlog(aT), while 
the Casimir energy is linear in T. However this linear term is independent 
of geometry, and thus, the energy difference between two balls of different 
radii decays exponentially. This is expected from the equipartition property 
of classical statistical mechanics - at high temperatures the average energy 
per mode, ksT, is independent of the frequency, and since there is a one to 
one correspondence between the eigenmodes of two spheres of different radii, 
the difference in the energy should vanish in the high temperature limit. The 



calculations were based on the method developed in [|14|. This method have 
since been used in other related problems at zero temperature in |T7, |TH], |T^. 
We expect that these cases now too admit an exact treatment, along similar 
lines as presented here. 
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